Análisis de Regresión con R

Comparación de modelos de regresión con R

Artículo científico

Datos

• Datos tomados del artículo referenciado previamente.
• Fuente de datos
• Nota: los análisis no arrojan exactamente los mismos resultados porque se tienen valores ausentes en algunas variables.

Lectura de datos

library(readxl)
library(tidyverse)
datos <- read_xlsx("../data/Data_Paper_Plos_One_Muscle.xlsx", skip = 3,
                   na = "N/A", n_max = 47) %>% 
  rename(RA_takeof_leg = RA...7,
         OB_takeof_leg = OB...8,
         PM_takeof_leg = PM...9,
         QL_takeof_leg = QL...10,
         ES_takeof_leg = ES...11,
         Gmax_takeof_leg = Gmax...12,
         Gmed_takeof_leg = Gmed...13,
         IL_takeof_leg = IL...14,
         RA_free_leg = RA...15,
         OB_free_leg = OB...16,
         PM_free_leg = PM...17,
         QL_free_leg = QL...18,
         ES_free_leg = ES...19,
         Gmax_free_leg = Gmax...20,
         Gmed_free_leg = Gmed...21,
         IL_free_leg = IL...22,
         id = ID,
         edad = `Age (years)`,
         altura_cm = `Height (cm)`,
         imc = `Body mass (kg)`,
         dist_salto_cm = `long jump distance (cm)`,
         sprint_100m_seconds = `100-m sprint time (s)`,
         grasa_subcut_cm2 = `Subcutaneous fat CSA (absolute value, cm2)`)
datos

Objetivos

• Replicar análisis estadísticos aplicados en el artículo científico de interés.
• Evidenciar la relación existente entre características anatómicas de atletas vs rendimiento en salto largo.
• Evaluar otros métodos de statistical learning y compararlos con los resultados obtenidos por los autores.

Resultados del paper

Correlaciones

• Aunque fueron numerosos los resultados obtenidos por los autores, para el objetivo de este documento se destacan los siguientes:
  • La relación entre el área transversal relativa (CSA) del recto abdominal (AR) del lado de la pierna de despegue y el mejor registro personal para el salto largo.
    • Correlación: 0.674
    • Valor p: 0.004 (estadísticamente significativo)
• Las correlaciones (con intervalo de confianza del 95%) se presentan en la siguiente tabla:

Gráfico de dispersión

Original

Réplica con R

library(ggplot2)
datos %>% 
  ggplot(data = ., aes(x = RA_takeof_leg, y = dist_salto_cm)) +
  geom_point(size = 3) +
  labs(x = expression('Relative CSA of RA takeoff leg side - cm'^"2"/'kg'^"2/3"),
       y = "Personal best record of long jump (cm)") +
  geom_smooth(method = "lm", se = FALSE, lty = 3, lwd = 1, color = "black") +
  theme_light()

Predichos vs Reales (paper)

Resultados adicionales con R

Distribuciones

Gráficos cuantil cuantil

Shapiro Wilk

Se comprueba la normalidad de las variables (\(\alpha = 0.05\)), bajo el siguiente juego de hipótesis:

\[H_0: X \sim N(\mu, \sigma^2)\\ H1: x \nsim N(\mu, \sigma^2)\]

Comparativos

• Se comparan registros de cross-sectional area (CSA) de la pierna de despeque (takeoff) vs la pierna libre (free). Las variables (músculos) a comparar son las siguientes:
  • RA: recto abdominal.
  • OB: oblicuos internos y externos.
  • PM: psoas mayor.
  • QL: cuadrado lumbar.
  • ES: erector spinae.
  • Gmax: gluteo mayor.
  • Gmed: gluteos medio y mínimo.
  • IL: iliaco

Matriz de correlaciones

• Se construye la matriz de correlaciones.
• La variable que presente mayor correlación lineal con la longitud del salto, será tenida en cuenta para estructurar el modelo de regresión lineal simple (RLS). Con las demás variables se construye el modelo de regresión lineal múltiple (RLM).
• Para la construcción del modelo de RLM se comprueba la multicolinealidad de las variables y se proponen tres alternativas:
  • Eliminación de variables por factor inflacionario de varianza.
  • Regresión por componentes principales.
  • Regresión regularizada (Ridge y Lasso)